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We systematically investigate the formation and stability of snowballs formed by femtosecond photo- 
ionization of small alkali clusters bound to helium nanodroplets. For all studied alkali species 
Ak=(Na, K, Rb, Cs) we observe the formation of snowballs Ak"'"HeN when multiply doping the 
droplets. Fragmentation of clusters Akn upon ionization appears to enhance snowball formation. 
In the case of Na and Cs we also detect snowballs Ak^HcN formed around Ak dimer ions. While 
the snowball progression for Na and K is limited to less than 11 helium atoms, the heavier atoms Rb 
and Cs feature wide distributions at least up to Ak~'"He4i. Characteristic steps in the mass spectra 
of Cs-doped helium droplets are found at positions consistent with predictions on the closure of the 
F^* shell of helium atoms around the Ak^ ion based on variational Monte Carlo simulations. 

PACS numbers: 32.80.Fb,36.40.Qv,82.50.Hp 



; I. INTRODUCTION 

Ion-induced electrostatic interactions are of fundamen- 
tal importance in various fields of physics, chemistry and 
biology |Jj, i2,]. Ions implanted in an ultracold helium 
(He) environment are of particular interest due to the 
influence of superfluidity on the mobility and on the He 
density distributions around the embedded ions. While 

I electron impurities reside in bubble states with low He 
density due to repulsive interaction, positive ions induce 
local enhancement of the He density by electrostriction. 
Since the He density in this compressed shell surround- 
ing the ion can exceed the density of solid He, such a 
structure is called a snowball. 

Using a phenomenological model by Atkins, the diam- 
eter of such a snowball in liquid He was estimated to 
about 6.5 A consisting of about 50 atoms in two He lay- 

I ers, independent of the specific element For alkali 
and alkaline earth metals a distinct dependence on the 

[ specific element was observed in ion mobility measure- 
ments 3. Compared to He'*^, mobilities were found to 
be lower for alkali and higher for alkaline earth atoms, in- 
dicating larger, respectively lower, numbers of He atoms 
attached to the ions. In addition, those experiments 
showed evidence that mobility increases with the mass 
of the dopant [3] . 

More recently, the stability of ion snowballs has been 
studied by means of abundance spectroscopy of snow- 

, balls formed upon ionization of metal atoms in He 
droplets [1, H, 0]. Snowballs with up to 150 He atoms 
attached to, e. g., Ag"'" and Mg+ ions have been observed. 
Element specific stability patterns reflecting closures of 
He solvation shells are directly observable as steps or 
kinks in the snowball abundance distributions. For ex- 
ample, closures of the first and second solvation shells 
around positive silver ions, Ag+, at Ni = 10 and N2=32 
He atoms, respectively, are clearly visible as steps in the 
mass spectrum after femtosecond (fs) photo-ionization of 
Ag clusters embedded in He droplets. A sensitive de- 
pendence of the abundance of snowball complexes on the 
ionizing laser intensity is observed. Upon strong field ion- 



ization followed by massive fragmentation of the metal 
kernel, snowballs with as many as 150 (20) He atoms at- 
tached to the ionic Ag"*" (Mg+) core appear in the mass 
spectra 0, Q. Moreover, the pump-probe dynamics of 
snowball formation was used as a diagnostic tool to inves- 
tigate the caging of metal clusters in large He droplets 0] . 

A number of theoretical techniques have been ap- 
plied to studying the solvation of positive ions in He 
droplets [1, 0|, [lo| . Mainly alkali and alkaline earth ions 
have been addressed so far since reliable Me """-He po- 
tentials are available for these species. Using variational 
Monte Carlo simulations, Reatto and coworkers find that 
all alkali and alkaline-earth cations form snowball struc- 
tures featuring shells of He atoms with high average den- 
sity [1, [§, [ll| . In addition to a modulated radial density 
profile around the impurity ions, snowballs are character- 
ized by angular correlations in the first He shell as well as 
by the degree of radial localization of He atoms and the 
ability of He atoms to move from the first to the second 
shell and vice versa. 

This solid-like order in the first shell turns out to be 
specific to the impurity ion. In contrast to the results of 
ion mobility measurements mentioned above, the num- 
ber of He atoms bound to the impurity ions in the first 
shell, Ni, increases with larger ion mass. E.g., for the 
alkali ions Na"*", K"*", and Cs"*", Rossi et al. have found 
Ni = 12, 15 and 17.5, respectively P. Most pronounced 
local order is found for the Na"*" ion which is explained by 
the deep Na+-He potential well and also by the match- 
ing of the minimum positions of the He-He and Na'''-He 
potentials. The calculated structure is compatible only 
with an icosahedron having 12 He atoms in the first shell. 
The Na+HcN system was also studied by Nakayama and 
Yamashita [l3|- They predicted a triple layered struc- 
ture with a rigid-body and solid-like first shell around 
the cation with 16 tightly bound He atoms in the first 
shell, considerably more than the 12 He atoms found in 
Ref. [|. 

The geometric structure of He shells in alkali-ion doped 
clusters was studied by quantum variational and diffusion 
Monte Carlo calculations [13, [H, [31 ■ By subsequently 
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adding He atoms to the ion and computing the result- 
ing cluster structure, marked differences in the AIc+Hcn 
geometry depending on N are found ■ The first solva- 
tion shell consists of equivalently arranged He atoms, i.e. 
their binding energy is determined by the two-body Ak+- 
He potential alone, independently from the other solvent 
atoms. Due to the screening of the charge by the He 
atoms in the first shell, the further atoms are bound only 
by the He-He interaction. During buildup of the first 
shell different regular structures are observed to emerge 
for Li+Hcg (octahedral), Li+Hcg and Li+Hcio, Na+Hcio 
and Na"'"Hei2 (icosahedral), K+Hcig and K"'"Hei2 (icosa- 
hedral). A sudden decrease of the binding energies per 
atom is found when adding an additional He atom to 
a snowball of regular structure. This is attributed to 
the newly-added atom pushing the already bound atoms 
away from the ion . The number of atoms in the first 
shell is found as Ni=10, 12 and 15 for Li+, Na+ and K+, 
respectively, in perfect agreement with the work cited 
above [8|. 

Similar calculations on He nanodroplets doped with 
cationic alkali dimers Akj'Hc N ( Ak = Li,Na, K) were 
performed by the same group jlSl . [iGj . Their results are 
characterized by the competition between the AkJ-He 
and the He-He potential. For Li^ the dimer-He interac- 
tion clearly dominates leading to a symmetric structure 
of all even-numbered clusters LijHe2N up to 2N=18. In 
the case of K J , on the other hand, the He atoms clearly 
arrange themselves in an asymmetric fashion, thus form- 
ing a helium mesh on one side of the dimer. Analysis 
of the cluster energetics yields values for the number of 
He atoms in the first solvation shell of Ni=6 for Li J and 
Na J , while in the K J case no sign of a shell closure was 
found. 

The study of the distribution of snowball sizes 
Ak^2HeN allows to check the predictions on the number 
of He atoms in the solvation shells around the ions. The 
present paper reports on an experimental study of the 
formation and stability of snowball complexes upon pho- 
toionization of alkali metal impurities doped onto helium 
nanodroplets using fs laser pulses. In contrast to other 
atoms or molecules, alkali atoms are weakly bound to the 
surface of He nanodroplets due to the extremely weak 
attractive interaction between Ak and He atoms p7j . 
Above-threshold laser ionization then leads to massive 
evaporation of He atoms as well as to the expulsion of 
the ionized dopant particle, which partly entrains He 
atoms to form a snowball structure. The experimental 
conditions for the formation of alkali helium snowball 
complexes are discussed with respect to the different be- 
haviour of monomer and dimer snowballs, Ak''"HeN and 
Ak^HcN- The experimental abundance spectra are com- 
pared to theoretical predictions, which agree with the 
experiment on the first He shell closure for Cs and Rb 
snowballs. These measurements extend earlier experi- 
ments on the formation and stability of alkali clusters em- 
bedded in helium nanodroplets using femtosecond photo 
ionization p^.[l9j. Furthermore, they provide important 



complementary information for the interpretation of fem- 
tosecond pump-probe experiments studying the dynam- 
ics of alkali- helium exciplex formation [20, [2l| and of the 
propagation of vibrational wave packets in alkali dimers 
formed on helium nanodroplets [22,, 23\ . 

II. EXPERIMENT 

In the experiment presented here, snowball complexes 
are studied by ionizing alkali (Ak) atoms and small clus- 
ters of various species attached to helium nanodroplets. 
The droplet machine employed for this study has already 
been described in detail elsewhere 0, [2^ . Briefiy, ultra- 
pure helium gas is expanded through an orifice 10 pm 
in diameter to form helium nanodroplets. The expan- 
sion conditions (pHo=60bar, nozzle temperature 13 K) 
are chosen within the subcritical expansion regime such 
that photo-ionization signals recorded both at the mass 
of the neat Ak+ ions as well as at the snowball masses 
are maximum. The corresponding average droplet size 
amounts to about 20,000 atoms. The source chamber 
is connected to the doping chamber by a Beam Dynam- 
ics skimmer, 400 pm in diameter. The He droplets are 
doped with Ak atoms by passing through a heated pickup 
cell before entering the detection volume of a commer- 
cial quadrupole mass spectrometer inside a third vacuum 
chamber. Here, the doped beam is ionized either by laser 
radiation or electron bombardment. 

We use a commercial modelocked fs Ti:Sa laser featur- 
ing a FWHM pulse length of 150 fs. The laser wavelength 
is selected for highest snowball yield relative to the neat 
ion Ak+ signal, thus resulting in different sets of wave- 
length and intensity for the different alkalis. The optimal 
wavelengths for Na, K, Rb and Cs are found to be 755, 
730, 800 and 860 nm, respectively. Using a lens with 
a focal length of / = 100 mm the peak laser intensity 
in the interaction volume was about 33, 23, 40 and 32 
GW/cm^, respectively. 

Upon photo-ionization, the ions Ak"*" desorb from 
the helium surface, partly entraining He atoms from 
the droplet to form snowball complexes. We find that 
the fragmentation of alkali clusters AkN upon photo- 
ionization into atomic ions Ak+ greatly enhances the for- 
mation of snowballs Ak+HcN- The clusters are formed by 
increasing the vapor pressure in the pickup cell, so that 
the helium droplets consecutively pick up more than a 
single atom. The temperature of the alkali oven was var- 
ied in order to maximize the snowball yield resulting in 
average loading of about 10 alkali atoms per droplet. 

III. RESULTS AND DISCUSSION 

We systematically study all Ak metals except Li, since 
the latter's small mass and isotopic composition seriously 
hamper the evaluation of mass spectrometric data of 
Ak+HcN snowball complexes. The snowball series 
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FIG. 1: Mass spectra displaying the Ak+HeN (Ak=Na, K) 
snowball progression on a logarithmic scale. Mass peaks re- 
sulting from contamination with other alkali metals or resid- 
ual water vapour are represented in gray. 

recorded in this experiment are depicted in figures[l]to|31 
Snowball formation can clearly be seen for all metals 
probed. The formation of small snowballs appears to 
be roughly equally efficient irrespective of the particu- 
lar Ak species considered: A comparison of the abun- 
dance ratios of the snowball Ak+He2 with respect to the 
ion Ak+ yields values between 0.5% and 1.0% for all 
shown mass spectra. The data obtained for Na, shown 
in Fig.[T]a), clearly reveal the presence of small snowballs, 
Na+Hci 3. Unfortunately the sample contains some per- 
centage of K, which — given a count rate for snowball 
formation of 0.5 % of the neat ion — is enough to super- 
impose a strong K mass peak onto the snowball series 
at the mass of Na"'"He4. All subsequent snowball peaks 
Na"''HeN>4 therefore coincide with the peaks of K snow- 
balls K+HeN- Nonetheless we are able to rule out the 
existence of snowballs larger than Na+Hcio under our 
experimental conditions. 

K behaves very similarly. The mass spectrum dis- 
played in Fig. [1] b) shows a snowball progression starting 
with the monomer and eventually disappearing in the 
dimer peak. Since the data displayed in figures [1] and [2] 
suffer from low resolution for technical reasons it is not 
possible to determine the exact end of the snowball se- 
ries. Nevertheless we can follow the peaks up to K+Hcs 
and can definitely rule out the presence of any snowballs 
larger than K+Heio. 

In contrast to K atoms, the heavy alkalis produce large 
snowballs with at least 41 helium atoms attached to the 
metal ion. Data obtained when doping Rb atoms are 
depicted in Fig. [51 The shaded peaks originate from Cs 
contamination of the sample and residual water vapor in 
the vacuum chamber. One can follow the snowball se- 
ries up to Rb"'"He4i; due to the isotopical broadening of 
the Rbjj" signal we can not unambigously assign further 
mass peaks to the monomer or the trimer progression 
Rb+HeN>4i or Rb;^HeN, respectively. Nonetheless, the 
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FIG. 2: Same as Fig. [TJ but for rubidium. The shaded area 
indicates the expected position of the first shell closure as 
derived from Q. Mass peaks resulting from contamination 
with cesium or residual water vapour are represented in gray. 

data clearly cover the mass range, in which the closure 
of the first He shell around the ionic core is expected 
according to variational Monte Carlo simulations [g] . Al- 
though Rossi et al. do not explicitly treat Rb"'"HeN snow- 
balls, their data show an increasing snowball size with 
Ak atom radius. Therefore, for Rb+ we expect a value 
Ni for the number of He atoms inside the first solvation 
shell between the value for K"*", Ni=15, and the one for 
Cs"*", Ni=17.5 He atoms. This range is indicated as a 
shaded area in Fig. O It is difhcult to see clear signs of 
a shell closure from the data presented in this figure, as 
the snowball distribution evolves smoothly through the 
region of interest. 

To unambiguously ascertain the position of shell clo- 
sure from our data, we plot the same data as 'differen- 
tial mass spectrum' in Fig. [3] a). In this representation 
we examine the abundance ratio of neighboring snow- 
ball mass intensities In/In-i, so this plot reveals pos- 
sible shell closures as dips: Supposing an exceptionally 
abundant snowball Ak"'"Heni, we expect In/In-i to drop 
significantly below unity for N=m-|-1. The data shown in 
this figure is the average over different measurements at 
various oven and nozzle temperatures to rule out a possi- 
ble influence of droplet or Ak cluster size on the snowball 
structure. The differential mass spectrum shows a sud- 
den drop for N=15, followed by a plateau. Compared to 
the simulations cited above the position of this drop al- 
most coincides with the calculated range of Ni=16-18.5. 

Fig. [4] features the mass spectrum obtained from Cs- 
doped He droplets. As in the case of Rb"*", we can follow 
the snowball distribution up to at least Cs"'"He4i. How- 
ever, there might well be more snowballs beyond the scan 
range employed for this measurement. In addition to the 
snowball series extending from the Cs^ monomer mass, 
we see clear evidence for snowballs formed around the Cs 
dimer ion, CsJHcn. The position of the closure of the 
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FIG. 3: 'Differential mass spectrum' of snowballs showing the 
relative abundances of neighbouring snowball masses AkHcN , 
AkHeN-i- As in figures [2] and |4l note the predicted position 
of the first shell closure. Some data points are missing due to 
an overlap of Cs, Cs2 peaks with the snowball masses. 

first shell, Ni=17.5 He atoms, as predicted in Ref. is 
also represented in the graph. In this mass range a step 
in the abundance spectrum of Cs^Hcn is clearly visible. 

To make sure the snowball distribution is not influ- 
enced by the ionization scheme employed in this mea- 
surement, additional mass spectra are taken by electron 
impact ionization. The resulting spectrum taken at an 
impact energy of 70 eV is shown in Fig. 2] as orange spec- 
trum, shifted by 1 amu towards higher masses. It shows a 
smooth decay of small snowball sizes Cs"''HeN<io in con- 
trast to the plateau found in the PI mass spectrum. More 
relevant, however, is the occurrence of the same step at 
Cs"'"Hei5...2o that was already found by PI. This shows 
that the step is indeed a feature intrinsic to the Cs+ ions 
independent of the ionization mechanism. 

The differential mass spectrum reveals another feature, 
that is easily overlooked in Fig. [D We note a sudden 
drop in intensity after Cs+Hei2, hinting at another rela- 
tively stable structure in addition to Cs+Heie. By means 
of Monte Carlo simulations of alkali-ion doped helium 
clusters, strong variations of the single-atom evaporative 
energy were found when adding an additional He atom 
to an ensemble of regularly arranged solvent atoms like 
the icosahedral K+Hei2 [12|. These variations are ex- 
plained by the large difference between the structures of 
the regular snowball Ak+HcN and the subsequent cluster 
Ak+HcN-i-i: the additional He atom pushes out the regu- 
larly bound solvent atoms to a position further away from 
the ion, thus decreasing their binding energy. Such a rel- 
atively weakly bound cluster is more susceptible to disso- 
ciation. Unfortunately there are no calculations available 
for Cs'''-doped helium clusters. However, it seems plau- 
sible to assume a similar behaviour for the insertion of a 
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FIG. 4: Same as Fig. [T] but for cesium. The Cs+Hcn pro- 
gression displays a distinct drop in snowball intensity in ac- 
cordance with the results in 0). Orange data were taken by 
EII and are shifted by 1 amu for comparison. 



He atom into Cs"'"Hei2, as seen in the mass spectrum in 
Fig.[3]b) for Cs+Hei3. 

It is worth mentioning that the spectra of Na and Cs 
also show dimer snowballs, Ak^HcN, while we are un- 
able to observe them for K and Rb. The Na progression 
Na^HcN can be traced to N=3 before disappearing in the 
mass peaks of KNa+HcN. However, since our study does 
not focus on snowballs around dimer ions Ak^HcN, they 
are strongly superimposed in the mass spectra by snow- 
balls formed around monomer ions Ak+ and impurities. 
Besides, the variations of the evaporation energy between 
neighbouring clusters are rather weak (ca. 25 cm^^ for 
Na^Hcg as compared to Na^Her) [TB]. The combination 
of these two facts might well account for the lack of a dis- 
tinct step in our snowball distribution. The same holds 
for our measurements of Cs2-doped snowballs, where the 
snowballs extend at least up to Cs^Hcg. Note that we 
cannot rule out the existence of larger snowballs, since 
they would exceed the mass range covered by our mea- 
surement. 

The formation of snowballs Ak^HeN was found to be 
governed by the relative strengths of the Ako'-He ver- 
sus the He-He interaction [l^. Since the Akg -He well 
is deepest for Li and substantially flattens out with in- 
creasing atomic number [24| . one could expect a mono- 
tonic change of snowball abundances when increasing the 
alkali dimer mass from NaJ to CsJ. However, this trend 
is probably masked by the effect of varying amounts of 
internal excitation carried by the dimer ions Ak^ of dif- 
ferent alkali species. Since a large fraction of the dimer 
ions AkJ originate from fragmentation of larger clusters, 
these AkJ may be produced in highly excited vibrational 
states. Thus, He atoms attached to Ak^ are evaporated 
to leave behind neat AkJ dimer ions, as previously ob- 
served by EII of rare-gas clusters in He droplets [l^ . Un- 
fortunately, it is beyond the scope of this investigation to 
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FIG. 5; Dependence of the yield of different snowballs on the 
amount of atoms doped into the droplets. Panel a): monomer 
snowballs Cs^HeN (N=0, 3, 6, 11), panel b): dimer snowballs 
CsjHeN (N=0, 3, 6). The gray dashed line indicates the atom 
density corresponding to most efficient doping of a single Cs 
atom. Note the logarithmic scale. 



assess the amount of energy carried by the individual al- 
kali dimer ions. 

The nonresonant nature of the ionization process is ex- 
pected to enhance fragmentation of the Cs clusters. This 
assumption is backed by PI spectra taken on snowball 
masses, which commonly show broad features; e.g., in 
the case of Cs^Hey we observe a broad maximum peak- 
ing around 755 nm with a FWHM of about 130 nm. The 
absence of the atomic Cs D line in the spectra indicates 
that even single-Cs snowballs Csj'^Hcn originate predom- 
inantly from Cs clusters. To assess the degree of frag- 
mentation occurring in our experiment, doping curves 
are recorded by varying the oven temperature and si- 
multaneously recording Cs^Hbn mass spectra. Typical 
examples of these curves are displayed in Fig. [51 with pan- 
els a) and b) showing the doping dependency of monomer 
snowballs Cs+Hcn and dimer snowballs CsJHcn, respec- 
tively. Within the plots the squares represent neat atoms 
or dimers Csf 21 while the circles, triangles and diamonds 
represent snowballs Cs^2-HeN- Taking into account the 
atom density for maximum doping of a single atom of 



2.1 • 10^^ cm~^ (indicated as a dashed line), we deduce 
that the vast majority of detected Csf 2 io'^s originate 
from fragments of clusters Csn formed on the droplets. 
Similar behaviour has been reported for the case of silver 
clusters [6]. 

It is interesting to note that the dependence of snow- 
ball abundance on the doping level differs from the one of 
neat Cs atoms: At Cs vapor pressures below 3-10^^ cm"'^ 
the curvature of the different doping curves increases with 
increasing N, indicating that bigger snowballs originate 
from bigger neat Cs clusters. However, further increas- 
ing the doping level above ca. 10^^ cm"-^ leads to falling 
snowball abundances before reaching the maximum in 
neat Cs intensity. This decline in the production of snow- 
balls with increasing alkali doping is related to the influ- 
ence of doping on the helium nanodroplet size: During 
the cluster growth each consecutive Cs atom dissipates 
binding energy into the droplet. This energy is released 
by evaporation of helium atoms, leading to a reduced 
number of helium atoms available for snowball forma- 
tion. 

This effect is more pronounced when photo-ionizing 
dimer snowballs. The abundance of dimer snowballs 
Cs^HeN (Fig- M b) passes its maximum at lower doping 
levels than the monomer snowball signal (Fig. [5] a). This 
is caused by the internal degrees of freedom of the ionized 
dimer Cs J . In addition to boiling off He atoms from the 
droplet upon ionization, the Cs2 bond can carry internal 
energy, which is dissipated by evaporative cooling of the 
snowball CsJHcn after it has left the droplet [25]. This 
leads to further reduction of the snowball abundance at 
high doping levels. 

Examining the maximum positions of the neat ions 
Cs+, CsJ, we find the monomer maximum to be shifted 
towards higher Cs doping than the dimer. This is in ac- 
cordance with observations of the fragmentation of alkali 
cluster ions Ak^ with Ak=Na, K. Fragmentation of large 
ionized clusters, Najlj and K^, predominantly produces 
monomer fragment ions Ak"*" as compared to dimer ions 

Ak+ [li,!!^. 

The doping curves recorded at monomer snowball 
masses Cs+Hcn (circles, triangles and diamonds in Fig. [5] 
a)) show a shoulder towards lower doping levels, which is 
very pronounced for small Cs^HeN<4, and whose width 
and intensity gradually fade away when observing larger 
snowballs up to CsHe2o . We take this as a sign of massive 
influx into these snowball masses from individual neat Cs 
cluster sizes. Based on the position of the peak relative 
to the position of maximal monomer doping, we roughly 
assign clusters around CS25 as the most probable precur- 
sor. 



IV. CONCLUSION 

In a systematic study of helium snowballs formed 
around alkali ions we have multiply loaded helium nano- 
droplets with all alkali metals except Li and created 
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snowballs by fs PI and by EII of alkali clusters. The data 
for Na and K allow only the study of rather small snow- 
balls of size up to Na+Hea and K+Heio- The results for 
the heavier alkali species permit comparison with varia- 
tional Monte Carlo simulations [Sj, ,12,] . The number of 
atoms observed in the first helium shell around Rb+ and 
Cs"*" ions is found to agree with theoretical predictions. 
The differential spectrum for Cs+ hints at an additional 
subshell closure for Cs"'"Hei2, possibly due to the swelling 
of the isocahedral helium shell around the Cs ion to ac- 
commodate the next He atom; however, for lack of theo- 
retical data this assumption cannot be verified. Snowball 
formation around alkali dimer ions is observed for NaJ 
and CsJ. 

Doping curves reveal strongly increased probability of 
formings snowballs around Cs+ and CsJ fragment ions 
coming from fragmenting large Cs clusters. The order 
of the maxima of doping curves measured for different 
snowball masses is interpreted in terms of massive evap- 
oration of helium upon cluster formation and ionization. 

In future experiments we plan to replace the fs multi- 
photon ionization scheme by single-photon ionization 
with nanosecond lasers. This will allow an exact deter- 
mination (and adjustment) of the energy deposited into 



the clusters and will thus drastically reduce fragmenta- 
tion. Furthermore, ionization just above the IP will rule 
out any contribution from higher excited ionic states. In 
that way the interpretation of our data will be greatly 
simplified. Furthermore, experiments with amplified fs 
pulses may induce much more extended snowball series, 
as observed with magnesium Fs pump-probe mea- 
surements may reveal the presence of caging-effects @ 
and shed some light onto the issue whether alkali clus- 
ters are surface bound or immersed into the He droplets. 
Furthermore, the theoretical study of snowball formation 
around dimer ions AkJ for the heavy alkalis Rb and Cs 
is desirable to interpret their element-specific appearance 
and stability patterns. Besides, modelling the dynamics 
of fragmentation of alkali clusters Akw on the droplets 
may add to the understanding of snowball formation. 
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